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Abstract 

It is well known that, for chaotic systems, the production of relevant entropy 
(Boltzmann-Gibbs) is always linear and the system has strong (exponential) sensi- 
tivity to initial conditions. In recent years, various numerical results indicate that 
basically the same type of behavior emerges at the edge of chaos if a specific gener- 
alization of the entropy and the exponential are used. In this work, we contribute to 
this scenario by numerically analysing some generalized nonextensive entropies and 
their related exponential definitions using z-logistic map family. We also corroborate 
our findings by testing them at accumulation points of different cycles. 
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1 Introduction 



In nonlinear dynamics, it is well known that the chaotic systems has an expo- 
nential sensitivity to initial conditions, characterized by the sensitivity func- 
tion (for one-dimensional case) 

X Ax(t) ,^ , 

iit) = lim — -H-> 1 



(where Ax{t) is the distance, in phase space, between two copies at time t) 
to diverge as ^{t) = exp(Ait), where Ai is the standard Lyapunov exponent 
[1]. On the other hand, for the marginal case, where Ai = 0, the form of the 
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sensitivity function could be a whole class of functions. One of the candidates 
is a power-law behavior, which can be characterized by an appropriate gen- 
eralization of exponentials, namely ^{t) = exp(At), where A is the generalized 
Lyapunov exponent. Therefore, A > and A < cases correspond to weak sen- 
sitivity and weak insensitivity to initial conditions respectively. This, in fact, 
constitutes a unified framework since the generalized exponentials include the 
standard one as a special case for an appropriate choice of related parameter. 

The other concept that we focus in this work is the entropy production. For 
a chaotic system, the Kolmogorov-Sinai (KS) entropy Ki is defined as the 
increase, per unit time, of the standard Boltzmann-Gibbs entropy and it is 
basically related to the standard Lyapunov exponents through the Pesin iden- 
tity, which states that Ki = Ai if Ai > and Ki = otherwise. Here, it is 
worth mentioning that the KS entropy is basically defined in terms of a single 
trajectory in phase space, using a symbolic representation of the regions of a 
partitioned phase space (sec [1]). However, it appears that, in most cases, this 
definition can be replaced by one based on an ensemble of initial conditions, 
which is the version we use herein. In the framework of this version, it has 
already been shown that the statistical definition of entropy production rate 
exhibits a close analogy to the production rate of thermal entropy and prac- 
tically coincides with the KS entropy in chaotic systems [2]. In recent years, 
there have been efforts on extending this picture to dynamical systems at their 
marginal points (like chaos threshold) by using a generalized entropic form, 
which allows us to define a generalized KS entropy K. From this definition, 
one can also conjecture a Pesin-like identity as K — X, which recovers both 
the chaotic and critical (i.e., chaos thresholds) cases since it includes standard 
Pesin identity as a special case [3] . 

As a whole, this unified framework has been numerically verified firstly for the 
logistic map [3,4] using the Tsallis entropy Sq = (l — Yli^ipfj / (q — 1) [5], 
which grows linearly for a special value of entropic index q, which is g^en — 0.24; 
whereas the asymptotic power-law sensitivity to initial conditions is character- 
ized by the generalized exponential exp(x) = expg(x) = [1 + (1 — q)x]^^^^~'^'^ 
with the same value of the entropic index. After these works on the logis- 
tic map, numerical evidences supporting this framework came also from the 
studies of other low-dimensional dynamical systems, such as the 2;-logistic map 
family [6,7], the Henon map [8] and the asymmetric logistic map family [9]. 
Besides these numerical investigations, analytical treatment of the subject is 
also available recently in a series of paper by Baldovin and Robledo [10]. Very 
recently, a similar analysis has been performed for the 2;-logistic map family, 
but this time, using ensemble-averaged initial conditions distributed uniformly 
over the entire available phase space [11,12]. The most important outcome of 
this analysis is another numerical verification of the coincidence of the entropic 
indices coming from the sensitivity function and entropy production rates (al- 
though with a different value ~ 0.36), which consequently broadens the 
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validity of the Pesin-like identity. 



Finally, here we should mention a recent work of Tonelli et al. [13], where 
they demonstrate that the above-mentioned framework is even more general 
by making use of a two-parameter family of logarithms [14] 

a -|- p 



where a characterizes the large (small) argument asymptotic behavior. 
Prom this wide class, they analysed four interesting one-parameter cases, 
namely ; 

(i) the Tsallis logarithm [5] : a = 1 — g and /3 = 

(ii) the Abe logarithm [15] : a = \ — q and /3 = q;/(1 -|- a) 

(iii) the Kaniadakis logarithm [16] : a — (3 — k — 1 — q 

(iv) the 7 logarithm : a — 1 — q and P — {1 — q)/2. 

Their analysis consists of studying the sensitivity function and the entropy 
increase rates for the logistic map at the edge of chaos. Obviously, for the 
corresponding entropy in each case, one needs to use 

S(t)^hMt)ii{-L.) (3) 



from where the entropy production rates in time can be calculated. Their 
numerical results clearly verified that, for the logistic map, the relavant value 
of the entropic index is g^g„ ~ 0.36 not only for the Tsallis case but also for 
the others as well, and moreover that the Pesin-like identity is also present for 
all cases (with different numerical values for each case). 

The aim of the present effort is to check the validity of the above-mentioned 
picture making use of the 2;-logistic map family. Moreover, the effect of different 
cycles on this vahdity is tested by analysing four distinct cycles. 



2 The model and the procedure 



The model system that we use in our analysis is the 2;-logistic map family 
defined as 

Xt+i ^ 1 - alxtl"" (4) 

where 2; > 1; < a < 2; \xt\ < 1; t = 0, 1,2, .... It is easily seen that z = 2 
case corresponds to the standard logistic map. 
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Table 1 

The critical Oc values of the z-logistic map family for all cycles used in this study. 



z 




a 


c 




cycle 2 


cycle 3 


cycle 4 


cycle 5 


1.75 


1.35506... 


1.74730... 


1.92764... 


1.60749... 


2 


1.40115... 


1.77981... 


1.94217... 


1.63101... 


2.5 


1.47054... 


1.82886... 


1.96144... 


1.66954... 


3 


1.52187... 


1.86299... 


1.97302... 


1.69944... 


4 


1.59490... 


1.90597... 


1.98524... 


1.74282... 



Firstly we study the sensitivity to initial conditions at the edge of chaos using 
the sensitivity function given in Eq.(l). Prom its definition, for the calcula- 
tion, we proceed with considering two initially very close points, which makes 
for example Aa;(0) = 10^^^ and then at each time step we numerically cal- 
culate the sensitivity function. To make an ensemble average, this procedure 
is repeated many times starting from different initial conditions all chosen 
randomly within the phase space and an average is taken over all values of 
log(^). The special value of a (consequently, g"e„) is found as the value for 
which we obtain a linear time dependence of (|log(^)^ (t). Finally, it is clear 
that the slope of this curve gives us the generalized Lyapunov exponent A. 
In our analysis we concentrate on various values of z for four different cycles 
at their accumulation points denoted by Cc- These critical values are given in 
Table 1. 



Secondly we study the entropy production rates of the entropies related to four 
selected logarithms, namely we use Eq.(3). The procedure that we implement 
for the entropy production is the same as the one introduced firstly in [2] 
for chaotic systems and then used in other works related to this context. It 
consists of dividing the phase space in W equal intervals and putting randomly 
initial points in one of them. Then, one should trace the spread of initial 
points within the phase space and calculate the entropy S{t\ from where we 
can obtain the entropy production per unit time as 

Sit) 

K = lim lim lim . (5) 



We then repeat this many times starting from randomly chosen different inter- 
vals of the phase space and average the entropy S(t) over these experiments. 
Here, the special value of a (consequently, is obviously found as the value 
for which we obtain a linear time dependence of (S) (t). In addition to this. 
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the value of K can be obtained from the slope of this linear curve. 

From the previous results [11,12], it is known for the Tsallis case that, for 
the z-logistic map family, the proper a values (namely, g^g„) coming from the 
sensitivity function and entropy production coincide each other and moreover, 
for all z values, for the proper a value the Pesin-like identity is preserved as 
K = X. Then, very recently, for the standard logistic map (namely, z = 2 case), 
Tonelli et al. [13] showed that the picture is wider that this by considering the 
other three cases mentioned above. Now, our aim is to check the validity of 
this framework for the z-logistic map and for various cycles. 



3 Numerical results and conclusions 

In our simulations, for the sensitivity function analysis, for each cycle and 
z value we consider two very close points so that Aa;(0) = 10~^^ and then 
numerically calculate the sensitivity function ^ from its definition. In order to 
take the logarithm of ^, we use four different cases mentioned previously. We 
repeat this for 4 x 10^ uniformly distributed initial conditions to make an en- 
semble average of ^log(^)^ (t). Its time evolution is given for a representative 
value of z for cycles 3, 4 and 5 in Fig.l(a),(c) and (e) respectively. For each 
case (namely, Tsallis, Abe, Kaniadakis and 7 cases), the special value of a is 
determined by testing various a values until we obtain a linear time depen- 
dence. From the slope of each curve, the corresponding generalized Lyapunov 
exponents A can be calculated. All of the obtained results for a and A are given 
in Table 2. For a given z value of each cycle, it is seen that the special value 
of a is the same (within error bars) for all four generalized cases. However, 
generically the generalized Lyapunov exponents are different in each case. 

Then, we study the entropy increase rates. We divide the phase space into 
W equal cells and put N initial conditions into a randomly chosen cell. We 
then let the dynamics evolve in time and calculate the generalized entropies. 
In order to make an ensemble average, we repeat this procedure many times 
starting from randomly chosen different cells (1^/2 cells), which reduces also 
the fluctuations that appear at the edge of chaos. In all our simulations we use 
W = 10^ and N = lOW. For each cntropic form, we use the special value of a 
calculated from the sensitivity function and we obtain the linear coefficient K 
from the slope of the linear entropy increase rates as it is seen in Fig.l(b),(d) 
and (f). As it is clear from Table 2 that, for all z values of each cycle and 
for all entropic forms, K values coincide with corresponding A values (within 
numerical errors). 

Finally, it is worth mentioning that generically for the increasing values of z 
(depending also on the cycle), the observed fluctuations in the time evolution 
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Fig. 1. Time evolution of the sensitivity function averaged over 4 x 10^ uniformly 
distributed initial conditions for (a) cycle 3 and z = 4, (c) cycle 4 and 2 = 4, 
(e) cycle 5 and z = 4. Entropy increase rates with W = 10^ and N = WW are 
given in (b),(d) and (f) for the same cycle and z values chosen for the sensitivity 
function. In each case an average over W/2 boxes are performed. 

of the sensitivity function become to be more stronger. This is even very much 
stronger in Kaniadakis case than other three cases. Due to these fluctuations, 
we did not give any value for z = 3,4 cases of cycle 2 in Table 2. A visualization 
of these fluctuations is given in Fig. 2. The reason for this is the following : If 
a given initial condition has weak insensitivity then ^ value is smaller than 1, 
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Fig. 2. (a) Time evolution of the sensitivity function averaged over 4 x 10^ uniformly 
distributed initial conditions for cycle 2 and z = 4. 

which makes the value of each generalized logarithm negative. But, whenever 
the ^ value is closer to 0, then logarithms start to go very large negative values 
and the most rapidly diverging case is the Kaniadakis one. As a result, these 
very large negative values dominate the sum for the average which yield large 
fluctuations. 

Summing up, in this work, we have numerically analysed the sensitivity to 
initial conditions and entropy increase rates of z-logistic map family at the edge 
of chaos for various cycles. We have discussed the same generalized logarithms 
and their related entropies studied in [13] for the standard logistic map and 
shown that the framework is similar to that of the standard logistic map; 
namely, (i) the proper a value, which gives a linear time dependence of the 
sensitivity function and a linear entropy production, is the same for all four 
entropic forms for a given z value of a given cycle, (ii) the results corroborate 
a generalized Pesin-like identity K = X for each entropic form with a different 
numerical value. 
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Table 2 

For the generalized logarithms given in the text, three important quantities are 
hsted for cycles 2, 3, 4 and 5 : (i) the appropriate a value which gives a linear time 

dependence of ^log(^)^ (t), (ii) the value of the generalized Lyapunov exponent 

A which is the slope of the above-mentioned linear time dependence and (iii) the 

value of entropy increase rates K which is calculated from the slope of the linearly 
increased entropies coming from each of the generalized logarithms. For many cases, 
the error bars are ±0.01 and in a few cases, they are found as ±0.02. 



cycle 


z 


Tsallis 


Abe 


Kaniadakis 


7 


a. 


K 


A 


a. 


K 


A 


a. 


K 


A 


a 


K 


A 


2 


2 


0.64 


0.26 


0.27 


0.65 


0.17 


0.18 


0.65 


0.14 


0.15 


0.65 


0.19 


0.19 


2.5 


0.66 


0.28 


0.28 


0.67 


0.19 


0.19 


0.69 


0.15 


0.18 


0.67 


0.20 


0.20 


3 


0.67 


0.28 


0.28 


0.67 


0.18 


0.18 








0.67 


0.19 


0.19 


4 


0.71 


0.31 


0.31 


0.72 


0.21 


0.21 








0.71 


0.21 


0.21 


3 


1.75 


0.08 


0.47 


0.48 


0.10 


0.37 


0.37 


0.10 


0.36 


0.35 


0.09 


0.42 


0.41 


2 


0.12 


0.48 


0.49 


0.14 


0.35 


0.36 


0.14 


0.33 


0.34 


0.13 


0.40 


0.40 


2.5 


0.18 


0.47 


0.48 


0.19 


0.30 


0.31 


0.19 


0.28 


0.29 


0.19 


0.37 


0.38 


3 


0.22 


0.44 


0.45 


0.23 


0.29 


0.29 


0.23 


0.26 


0.26 


0.23 


0.35 


0.35 


4 


0.28 


0.41 


0.41 


0.29 


0.2G 


0.2G 


0.29 


0.23 


0.23 


0.28 


0.28 


0.28 


4 


1.75 


0.01 


0.58 


0.59 


0.03 


0.56 


0.57 


0.03 


0.56 


0.56 


0.02 


0.58 


0.60 


2 


0.02 


0.57 


0.57 


0.05 


0.53 


0.55 


0.05 


0.53 


0.54 


0.03 


0.55 


0.56 


2.5 


0.05 


0.56 


0.58 


0.07 


0.47 


0.47 


0.07 


0.46 


0.45 


0.06 


0.51 


0.52 


3 


0.07 


0.52 


0.53 


0.09 


0.41 


0.41 


0.09 


0.40 


0.39 


0.09 


0.48 


0.50 


4 


0.11 


0.46 


0.46 


0.13 


0.34 


0.34 


0.13 


0.32 


0.32 


0.13 


0.40 


0.42 


5 


1.75 


0.04 


0.40 


0.41 


0.06 


0.35 


0.35 


0.06 


0.34 


0.34 


0.05 


0.37 


0.38 


2 


0.07 


0.39 


0.42 


0.09 


0.32 


0.32 


0.09 


0.31 


0.31 


0.08 


0.35 


0.36 


2.5 


0.12 


0.37 


0.38 


0.14 


0.28 


0.28 


0.14 


0.26 


0.26 


0.13 


0.31 


0.31 


3 


0.16 


0.35 


0.35 


0.18 


0.25 


0.25 


0.18 


0.23 


0.24 


0.17 


0.28 


0.28 


4 


0.22 


0.30 


0.30 


0.23 


0.19 


0.19 


0.23 


0.18 


0.18 


0.23 


0.23 


0.23 
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